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Abstract 

Evidence for the hadronic annihilation decay mode B + — > Df(p is found with greater 
than 3d significance. The branching fraction and CP asymmetry are measured to be 

B{B + -> D+(f>) = (l.87to%t (stat) ± 0.19 (syst) ± 0.32 (norm)) x 10~ 6 , 
A C p(B + -> D+(f>) = -0.01 ±0.41 (stat) ±0.03 (syst). 

The last uncertainty on B(B + — > Df(f>) is from the branching fractions of the 
B + — > DfD° normalization mode and intermediate resonance decays. Upper limits 
are also set for the branching fractions of the related decay modes B^ — > D^K*°, 

B (c) -> D (s)K*° and B+ -> D+(f>, including the result B{B + D+K*°) < 1.8 x 10~ 6 
at the 90% credibility level. 
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1 Introduction 



The decaysQfi 4 " D+4>, D + K*°, D+K*° occur in the Standard Model (SM) via anni- 
hilation of the quarks forming the B + meson into a virtual W + boson (Fig. [T). There 
is currently strong interest in annihilation- type decays of B + mesons due, in part, to 
the roughly 2a deviation above the SM prediction observed in the branching fraction 
of B + — > r + u [lj|2]- Annihilation diagrams of B + mesons are highly suppressed in the 
SM; no hadronic annihilation-type decays of the B + meson have been observed to-date. 
Branching fraction predictions (neglecting rescattering) for B + —t Df(p and B + — >■ D + K*° 
are (1 — 7) x 10~ 7 in the SM [3]-[6], where the precision of the calculations is limited 
by hadronic uncertainties. The branching fraction for the B + — > D+K*° decay mode 
is expected to be about 20 times smaller due to the CKM quark-mixing matrix ele- 
ments involved. The current upper limits on the branching fractions of these decay 
modes are B{B + ->• D+cj)) < 1.9 x 1(T 6 |7|, B(B + -»■ D+K*°) < 3.0 x 10~ 6 g and 
B(B + ^ D+K*°) < 4.0 x 10" 4 [9], all at the 90% confidence level. 

Contributions from physics beyond the SM (BSM) could greatly enhance these branch- 
ing fractions and/or produce a large CP asymmetry [3]. For example, a charged Higgs 
(H + ) boson mediates the annihilation process. Interference between the W + and H + 
amplitudes could result in a CP asymmetry if the two amplitudes are of comparable size 
and have both strong and weak phase differences different from zero. An H + contribution 
to the amplitude could also significantly increase the branching fraction. 

In this paper, first evidence for the decay mode > D+cf) is presented using 1.0 fb _1 
of data collected by LHCb in 2011 from pp collisions at a center-of-mass energy of 7 TeV. 
The branching fraction and CP asymmetry are measured. Limits are set on the branching 
fraction of the decay modes B + — > D + K*° and B + — > DfK*°, along with the highly 
suppressed decay modes B + — > D + K*° and B + — > DfK*°. Limits are also set on the 
product of the production rate and branching fraction for B+ decays to the final states 
D+(j>, D+ s) K*° and D^K*°. 

2 The LHCb experiment 

The LHCb detector (To] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes placed downstream. The combined tracking 
system has a momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 
lOOGeV/c, and an impact parameter resolution of 20 (a.m for tracks with high transverse 
momentum (pr)- Discrimination between different types of charged particles is provided 

throughout this paper, charge conjugation is implied. Furthermore, K*° and (f> denote the if* (892) 
and 0(1020) resonances, respectively. 
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by two ring- imaging Cherenkov detectors 11 . Photon, electron and hadron candidates are 



identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
muon system composed of alternating layers of iron and multiwire proportional chambers. 

The LHCb trigger 12 consists of a hardware stage, based on information from the 
calorimeter and muon systems, followed by a software stage which applies a partial event 
reconstruction (only tracks with px > 0.5GeV/c are used). The software stage of the 
LHCb trigger builds two-, three- and four-track partial 6-hadron candidates that are 
required to be significantly displaced from the primary interaction and have a large sum 
of pr i n their tracks. At least one of the tracks used to form the trigger candidate must 
have > 1.7GeV/c and impact parameter x 2 with respect to the primary interaction 
Xip > 16. The Xip i s defined as the difference between the x 2 °f the primary interaction 
vertex reconstructed with and without the considered track. A boosted decision tree 
(BDT) fl3 14 is used to distinguish between trigger candidates originating from b- hadron 



decays and those that originate from prompt c-hadrons or combinatorial background. The 
BDT provides a pure sample of bb events for offline analysis. 

For the simulation, pp collisions are generated using Pythia 6.4 [15] with a specific 
LHCb configuration 

in which final state radiation is generated using Photos 
generated particles with the detector and its response are implemented using the Geant4 
toolkit 



Decays of hadronic particles are described by EvtGen [17] 



18 . The interaction of the 



19 as described in Ref. |20|. 



3 Event selection 

Candidates of the decays searched for are formed from tracks that are required to have 
Pt > 0.1 GeV/c, Xj P > 4 and p > 1 GeV/c. For the and K*° decay products the mo- 
mentum requirement is increased to p > 2 GeV/c. The Df — > K + K~n + , D + — > K~7t + tt + , 
cf) — > K + K~ and K*° — > K + tt~ candidates are required to have invariant masses within 
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25, 25, 20 and 50MeV/c 2 of their respective world-average (PDG) values j2TJ. The decay 
chain is fit constraining the Df, candidate mass to its PDG value. The Df, vertex is 

to (s) (s) 

required to be downstream of the B + vertex and the p-value formed from xfp + Xvertcx °f 
the B + candidate is required to be greater than 0.1%. Backgrounds from charmless decays 
are suppressed by requiring significant separation between the Df, and B + decay vertices. 

Cross-feed between D + and Df candidates can occur if one of the child tracks is 
misidentified. If a Df — > K + K~7r + candidate can also form a D + — > K~tt + tt + candi- 
date that falls within 25 MeV/c 2 of the PDG D + mass, then it is rejected unless either 
\ m KK — m DG | < 10 MeV/c 2 or the ambiguous child track satisfies a stringent kaon particle 
identification (PID) requirement. This reduces the D + — > Df cross-feed by a factor of about 
200 at the expense of only 4% of the signal. For decay modes that contain a D + meson, a 
D + — > K~7t + tt + candidate that can also form a Df — > K~K + n + candidate whose mass is 
within 25MeV/c 2 of the PDG Df mass is rejected if either \m KK — m^ DG | < 10MeV/c 2 
or the ambiguous child track fails a stringent pion PID requirement. For all modes, 
A+ — > Df, cross-feed (from the A+ — > pK~n + decay mode) is suppressed using similar 
requirements. 

When a pseudoscalar particle decays into a pseudoscalar and a vector, V, the spin 
of the vector particle (in this case a or K*°) must be orthogonal to its momentum to 
conserve angular momentum; i.e., the vector particle must be longitudinally polarized. For 
a longitudinally-polarized (K*°) decaying into the K + K~(K + ir~) final state, the angular 
distribution of the K + meson in the V rest frame is proportional to cos 2 9k, where 9 k is 
the angle between the momenta of the K + and B + in the V rest frame. The requirement 
| cos 6^ | > 0.4, which is 93% efficient on signal and rejects about 40% of the background, 
is applied in this analysis. 

Four BDTs that identify Df K + K'tt + , D + -»■ K~<k + <k + , K+K~ and 
K*° — > K + tt~ candidates originating from 6-hadron decays are used to suppress the 
backgrounds. The BDTs are trained using large clean Dfx, and K*° samples obtained 

from B ^ — > D^n~ , B® — > J/ip(fi and B° — > J/ipK*° data, respectively, where the back- 
grounds are subtracted using the sPlot technique (22). Background samples for the training 
are taken from the Df,, <p and K*° sidebands in the same data samples. The BDTs 
take advantage of the kinematic similarity of all 6-hadron decays and avoid using any 
topology-dependent information. The BDTs use kinematic, track quality, vertex and PID 
information to obtain a high level of background suppression. In total, 23 properties per 
child track and five properties from the parent Df^, <fi or K*° meson are used in each BDT. 
The boosting method used is known as bagging |23] , which produces BDT response values 
in the unit interval. 



*o 



A requirement is made on the product of the BDT responses of the Df, and or K 
candidates. Tests on several B9, — > DD' decay modes show that this provides the best 



performance 24 . The efficiencies of these cuts are obtained using large B® s , — > Df,7r , 
B° s J/0 and B° — > J/ip K*° data samples that are not used in the BDT training. The 
efficiency calculation takes into account the kinematic differences between the signal and 
training decay modes using additional input from simulated data. Correlations between 
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Table 1: Summary of fit regions for B + — >■ D^~<f). About 89% of the signal is expected to 
be in region A. 





\m KK 


-77ty| (MeV/c 2 ) 


cos 6 k 


< 20 


(20,40) 


> 0.4 


A 


B 


< 0.4 


C 


D 



the properties of the and or K*° mesons in a given B + candidate are also accounted 
for. 

The optimal BDT requirements are chosen such that the signal significance is maximized 
for the central value of the available SM branching fraction predictions. The signal efficiency 
of the optimal BDT requirement is 51%, 69% and 51% for B + D+<f>, B + -»■ D + K*° and 
B + — > DfK*° decay modes, respectively. The final sample contains no events with multiple 
candidates. Finally, no consideration is given to contributions where the K + K~(K + n~) 
is in an S'-wave state or from the tails of higher (f)(K*°) resonances. Such contributions 
are neglected as they are expected to be much smaller than the statistical uncertainties. 

4 Branching fraction for the B + — > D^4> decay 

The B + — > yield is determined by performing an unbinned maximum likelihood fit to 
the invariant mass spectra of B + candidates. Candidates failing the cos8k and/or tukk 
selection criteria that are within 40 MeV/ c 2 of m™ G are used in the fit to help constrain 
the background probability density function (PDF). The data set is comprised of the four 
subsamples given in Table [1} They are fit simultaneously to a PDF with the following 
components: 

• B + — > Df<p: A Gaussian function whose parameters are taken from simulated data 
and fixed in the fit is used for the signal shape. The fraction of signal events in 
each of the subsamples is also fixed from simulation to be as follows: (A) 89%; (B) 
4%; (C) 7% and (D) no signal expected. Thus, almost all signal events are expected 
to be found in region A, while region D should contain only background. A 5% 
systematic uncertainty is assigned to the branching fraction determination due to 
the shape of the signal PDF. This value is obtained by considering the effect on the 
branching fraction for many variations of the signal PDFs for B + — > D^ifi and the 
normalization decay mode. 

• B + — > D* + (p: The in this decay mode does not need to be longitudinally polarized. 
When the photon from the D* + decay is not reconstructed, the polarization affects 
both the invariant mass distribution and the fraction of events in each of the 
subsamples. Studies using a wide range of polarization fractions, with shapes taken 
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from simulation, show that the uncertainties in this PDF have a negligible impact 
on the signal yield. 

B° s K~ K*°: These decay modes, which arise as backgrounds to B + — > 

when the pion from the K*° decay is not reconstructed, have not yet been observed; 
however, they are expected to have similar branching fractions to the decay modes 
B° D^ + K-K*°. The ratio B(B° S D*+K~ K*°)/B(B° S D+K-K*°)js fixed 
to be the same as the value of B(B° D*+K-K*°)/B(B° D+K~K*°) 
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The 

fraction of events in each subsample is constrained by simulation. Removing these 
constraints results in a 1% change in the signal yield. 

• Combinatorial background: An exponential shape is used for this component. The 
exponent is fixed to be the same in all four subsamples. This component is assumed 
to be uniformly distributed in cos 6^. Removing these constraints produces shifts in 
the signal yield of up to 5%; thus, a 5% systematic uncertainty is assigned to the 
branching fraction measurement. 

To summarize, the parameters allowed to vary in the fit are the signal yield, the yield and 
longitudinal polarization fraction of B + — > D* + <p, the yield of B° s — > Di*^ + K~ K*° in each 
subsample, the combinatorial background yield in each subsample and the combinatorial 
exponent. 

Figure [2] shows the B + candidate invariant mass spectra for each of the four subsamples, 
along with the various components of the PDF. The signal yield is found to be 6.7+2 % where 
the confidence interval includes all values of the signal yield for which log (£ m ax/£) < 0.5. 
The statistical significance of the signal is found using Wilks Theorem [26] to be 3.6<j. A 
simulation study consisting of an ensemble of 10 5 data sets confirms the significance and 
also the accuracy of the coverage to within a few percent. All of the variations in the 
PDFs discussed above result in significances above 3a; thus, evidence for B + — > Df(f> is 
found at greater than 3a significance including systematics. 

The B + — > Df(p branching fraction is normalized to B(B + — >■ DfD°). The selection 
for the normalization mode, which is similar to that used here for B + — > Dt<b, is described 



in detail in Ref. 24 . The ratio of the efficiency of the product of the geometric, trigger, 
reconstruction and selection (excluding the charmless background suppression and BDT) 
requirements of the signal mode to the normalization mode is found from simulation to be 
0.93 ± 0.05. The ratio of BDT efficiencies, which include all usage of PID information, 
is determined from data (see Sect. [3j) to be 0.52 ± 0.02. The large branching fraction of 
the normalization mode permits using a BDT requirement that is nearly 100% efficient. 
For the charmless background suppression requirement, the efficiency ratio is determined 
from simulation to be 1.15 ± 0.01. The difference is mostly due to the fact that the 
normalization mode has two charmed mesons, while the signal mode only has one. The 
branching fraction is measured as 

B(R+^ D + rh) e(B+^DtD°)B(D°^K-n+) N(B+^ Df<P) + - 
B(B ^ Ds<P) ~ e{B+^Di<P) m -+ K+K=) N(B+^ D^W) ^ ~+ ° s ° ] 
= (l.87l^ (stat) ± 0.19 (syst) ± 0.32 (norm)) x 10~ 6 , 
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Figure 2: Fit results for B + — >■ The fit regions, as given in Table [T| are labelled on 

the panels. The PDF components are as given in the legend. 



where e denotes efficiency. The normalization uncertainty includes contributions from 
B(B + ^ D+D°) = (1.0±0.17)%, B(D° K~Tr + ) = (3.88±0.05)% and B{(j) K+K~) = 
(48.9 ± 0.5)% [2l]. The systematic uncertainties are summarized in Table [2] The value 
obtained for B(B + — > D+<fi) is consistent with the SM calculations given the large uncer- 
tainties on both the theoretical and experimental values. 
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Table 2: Systematic uncertainties contributing to B(B + — > Df<f))/B(B + DfD°). 



Source 


Uncertainty (%) 


Selection 


7 


Signal PDF 


5 


Background PDF 


5 


Normalization 


17 



5 Branching fractions for the decays B + — > D^K*° 
and B + 23+ K*° 

The SM predicts the branching fraction ratios B(B + D + K*°)/B(B + D+(p) ~ 
1 and B(B+ -»■ D+K*°)/B(B + £>+0) ~ |Kd/K s | 2 N- The partially reconstructed 



backgrounds are expected to be much larger in these channels compared to B + — > 
mainly due to the large K*° mass window. Producing an exhaustive list of decay modes 
that contribute to each of these backgrounds is not feasible; thus, reliable PDFs for the 
backgrounds are not available. Instead, data in the sidebands around the signal region are 
used to estimate the expected background yield in the signal region. The signal region 
is chosen to be ±2cr around the B + mass, where a = 13.8MeV/c 2 is determined from 
simulation. 

Our prior knowledge about the background can be stated as the following three 
assumptions: (1) the slope is negative, which will be true provided &-baryon background 
contributions are not too large; (2) it does not peak or form a shoulder^] and (3) the 
background yield is non-negative. These background properties are assumed to hold 
throughout the signal and sideband regions. To convert these assumptions into background 
expectations, ensembles of background-only data sets are generated using the observed 
data in the sidebands and assuming Poisson distributed yields. For each simulated data 
set, all interpolations into the signal region that satisfy our prior assumptions are assigned 
equal probability. These probabilities are summed over all data sets to produce background 
yield PDFs, all of which are well described by Gaussian lineshapes (truncated at zero) 
with the parameters /ibkgd and c bkgd given in Table [3j The B + candidate invariant mass 
distributions, along with the background expectations, are shown in Fig. |3j The results 
of spline interpolation using data in the sideband bins, along with the 68% confidence 
intervals obtained by propagating the Poisson uncertainties in the sidebands to the splines, 
are shown for comparison. As expected, the spline interpolation results, which involve a 
stronger set of assumptions, have less statistical uncertainty. 



A Bayesian approach 27 is used to set the upper limits. Poisson distributions are 



assumed for the observed candidate counts and uniform, non-negative prior PDFs for the 

2 No evidence of peaking backgrounds is found in either the Df, or K*° sidebands. If peaking 
backgrounds do make significant contributions, then the limits set in this paper are conservative. 
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Figure 3: Invariant mass distributions for (a) B + ->■ D + K*°, (b) B + ->■ D + K*°, 
(c) fi+ -»■ D+iT* and (d) fi + -»■ DfK*°. The bins are each 4cx wide, where 
a = 13.8 MeV/c 2 is the expected width of the signal peaks (the middle bin is centred at the 
expected B + mass). The shaded regions are the /ibkgd ± o"bk g d intervals (see Table [3]) used 
for the limit calculations; they are taken from the truncated- Gaussian priors as discussed 
in the text. Spline interpolation results (solid blue line and hashed blue areas) are shown 
for comparison. 



signal branching fractions. The systematic uncertainties in the efficiency and B + — > DfD° 
normalization are encoded in log-normal priors, while the background prior PDFs are the 
truncated Gaussian lineshapes discussed above. The posterior PDF, p(B\n ohs ), where n ohs 
is the number of candidates observed in the signal region, is computed by integrating over 
the background, efficiency and normalization. The 90% credibility level (CL) upper limit, 
Bqq, is the value of the branching fraction for which J® 90 p(B\n \ iS )dB = 0.9 J °° p(B\n ^ s )dB. 
The upper limits are given in Table [3j The limit on B + — > D + K*° is 1.7 times lower 
than any previous limit, while the B + — > DfK*° limit is 91 times lower. For the highly 
suppressed decay modes B + — > D + K*° and B + — y D+K*° these are the first limits to be 
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Table 3: Upper limits on B(B ± — > D,sK*°), where n Q b s is the number of events observed 
in each of the signal regions, while /ibkgd and <7bkgd are the Gaussian parameters used in 
the background prior PDFs. 



Decay 


n b s 


/^bkgd 


Obkgd 


Upper Limit at 90% CL 


B + ^D + K*° 


8 


2.2 


3.4 


1.8 x 10~ 6 


B + ^D + K*° 


8 


7.1 


3.6 


1.4 x 10~ 6 


B + ->• D+K*° 


19 


20.0 


4.2 


3.5 x 10~ 6 


B + D+K*° 


16 


14.8 


5.6 


4.4 x 10~ 6 



set. 

The posterior PDF for the B + — > D + K*° decay excludes the no-signal hypothesis at the 
89% CL and gives a branching fraction measurement of B(B + — > D + K*°) = (0.8 IqJ) x 10~ 6 , 
where the uncertainty includes statistics and systematics. This result is consistent with 
both the SM expectation and, within the large uncertainties, with the value obtained 
above for B(B + — > -D+0). If processes beyond the SM are producing an enhancement 
in B(B + — >■ Df(p), then a similar effect would also be expected in B + — > D + K*°. While 
an enhancement cannot be ruled out by the data, the combined B(B + — > -D+ 0) and 
B{B + — > D + K*°) result is consistent with the SM interpretation. 

6 Limits on branching fractions of decay modes 

Annihilation amplitudes are expected to be much larger for f?+ decays due to the large 
ratio of iV^/Ki&l- In addition, the B+ — > -D+0, D + K*°, D+K*° decay modes can also 
proceed via penguin-type diagrams. However, due to the fact that B+ mesons are produced 
much more rarely than B + mesons in 7 TeV pp collisions (the ratio of B£ to B + mesons 
produced is denoted by j c j f u ), no signal events are expected to be observed in any of 
these B^t channels. The Bayesian approach is again used to set the limits. A different 
choice is made here for the background prior PDFs because the background levels are so 
low. The background prior PDFs are now taken to be Poisson distributions, where the 
observed background counts are obtained using regions of equal size to the signal regions 
in the high-mass sidebands. Only the high-mass sidebands are used to avoid possible 
contamination from partially reconstructed B^ backgrounds. In none of the decay modes 
is more than a single candidate seen across the combined signal and background regions. 
The limits obtained, which are set on the product of f c /f u and the branching fractions 
(see Table [4]) , are four orders of magnitude better than any previous limit set for a B+ 
decay mode that does not contain charmonium. As expected given the small numbers of 
candidates observed, the limits have some dependence on the choice made for the signal 
prior PDF. As a cross check, the limits were also computed using various frequentist 
methods. The largest difference found is 20%. 
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Table 4: Upper limits on f c /f u 4 B(B C — > X), where n Q b s and n bkgd are the number of 
events observed in the signal and background (sideband) regions, respectively. 



Decay 


^obs 


^bkgd 


Upper Limit at 90% CL 


B+ D+<j) 








0.8 x 10~ 6 


B+ ->■ D + K*° 


1 





0.5 x 10~ 6 


B+ -»■ D+i?* 








0.4 x 10~ 6 


£+ ->■ DfK*° 








0.7 x 10~ 6 


B+ -»■ D+K*° 


1 





1.1 x 10~ 6 



7 CP asymmetry for the decay B + — >■ D^4> 

To measure the CP asymmetry, Aqp, in B + — > D+cfr, only candidates in region (a) and 
in a ±2<r window (±26.4 MeV/c 2 ) around the B + mass are considered. The number of 
B + candidates is n + = 3, while the number of B~ candidates is n_ = 3. The integral 
of the background PDF from the fit described in detail in Sect. [4] in the signal region is 
n bkgd — 0.75 (the background is assumed to be charge symmetric). The observed charge 
asymmetry is Abs = — n +)/( n - +n + — nbkgd) = 0.00 ±0.41, where the 68% confidence 
interval is obtained using the Feldman-Cousins method |28|. 

To obtain Acp, the production, ^4 pro d, reconstruction, v4 reco , and selection, „4. sc i, 
asymmetries must also be accounted for. The D+4> final state is charge symmetric except 
for the pion from the Z?+ decay. The observed charge asymmetry in the decay modes 
B + —7- J/tpK + and B + — > D°tt + , along with the interaction asymmetry of charged 
kaons (29] and the pion-detection asymmetry (30] in LHCb are used to obtain the estimate 
^Iprod ± A-eco = (— 1 ± 1)%- The large B° s — > Dfir~ sample used to determine the BDT 
efficiency is employed to estimate the selection charge asymmetry yielding ^4 se i = (2 ± 3)%, 
where the precision is limited by the sample size. Finally, the CP asymmetry is found to 
be 

A CP (B + ^ D+cj)) = A ohs - A prod - Acco - Ad = -0.01 ± 0.41 (stat) ± 0.03 (syst), 
which is consistent with the SM expectation of no observable CP violation. 

8 Summary 

The decay mode B + — > D+<fi is seen with greater than 3a significance. This is the first 
evidence found for a hadronic annihilation-type decay of a B + meson. The branching 
fraction and CP asymmetry for B + — > Df<j) are consistent with the SM predictions. 
Limits have also been set for the branching fractions of the decay modes Bt, — > Dt,K*°, 

Bts -> D+K*° and £?+ -)■ D+cfi. These limits are the best set to-date. 
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